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Nucleon spin fluctuations in a dense medium reduce the "naive" values of weak interaction rates 
(neutrino opacities, neutrino emissivities) . We extend previous studies of this effect to the degenerate 
case which is appropriate for neutron stars a few ten seconds after formation. If neutron-neutron 
interactions by a one-pion exchange potential are the dominant cause of neutron spin fluctuations, 
a perturbative calculation of weak interaction rates is justified for T < 3m/(47ra^) » 1 MeV, where 
m is the neutron mass and a n » 15 the pion fine-structure constant. At higher temperatures, the 
application of Landau's theory of Fermi liquids is no longer justified, i.e. the neutrons cannot be 
viewed as simple quasiparticles in any obvious sense. 

PACS numbers: 97.60.Jd, 13.15.+g 14.60.Lm, 95.30.Cq 



I. INTRODUCTION 



In a dense nuclear medium the effective neutrino in- 
teraction rates are modified by the presence of nucleon- 
nucleon interactions. While the importance of spatial 
spin-spin correlations has been recognized for a long time, 
it had been overlooked that the interaction-induced tem- 
poral fluctuations of the spin of a single nucleon can be 
a more important effect. It reduces the naive neutrino 
opacities and neutrino emissivities of nuclear matter be- 
low their naive values jl]-§] . These studies focussed on a 
classical nucleon plasma, i.e. on nonrelativistic and non- 
degenerate conditions which are thought to obtain in the 
core of a supernove for the first few seconds after col- 
lapse. It was found that the spin-fluctuation rate in this 
environment is so large that it is not possible to calcu- 
late weak interaction rates by a perturbative expansion 
in terms of the nucleon-nucleon interaction potential. 

We presently study the same effect for a degenerate 
medium in order to derive a perturbative expression for 
the cross-section reduction by nucleon spin fluctuations, 
and in order to understand the physical conditions of 
temperature and density where a "naive" calculation of 
weak interaction rates may be justified. Many attempts 
have been made to calculate neutrino opacities and emis- 
sivities for the physical conditions pertaining to neutron 
stars because of the obvious importance of these quanti- 
ties for a theoretical understanding of neutron-star cool- 
ing §. While many of these works are dedicated to calcu- 
lating the impact of spatial correlation effects on neutrino 
interaction rates, none of them appears to have addressed 
the important issue of nucleon spin autocorrelations. 

One may take a somewhat different perspective on the 
same problem if one notes that these calculations are 
based on Landau's theory of interacting Fermi liquids 



where a "nucleon" is a quasiparticle excitation of the 
medium |6). This picture is justified only if the quasipar- 
ticles near the Fermi surface do not interact too strongly, 
i.e. r _1 <C T, where r is a typical time between collisions 
and T is the temperature of the medium. Landau's condi- 
tion is based on the observation that at T = the Fermi- 
Dirac distribution is a step function which, at nonzero 
temperature, is smeared out over an energy range of ap- 
proximate width T. Collisions, on the other hand, in- 
troduce an energy uncertainty of order r _1 which clearly 
should be much smaller than T in order for the Fermi- 
Dirac distribution to make any sense. When Landau's 
condition is violated it is not possible to speak of quasi- 
particles which obey Fermi statistics. Degeneracy effects 
ensure that the time between collisions becomes large at 
low temperatures, so there is no significant restriction in 
the T — > limit. For hot neutron-star matter, however, 
it is not a priori obvious that Landau's condition is sat- 
isfied. We were unable to locate any discussion of this 
problem in the entire literature pertaining to weak in- 
teraction rates in neutron stars [g. Therefore, it is not 
frivolous to raise the question of how cold the medium 
in a neutron star has to become before a Fermi-liquid 
treatment becomes possible. 

As previously argued [|]-[| , the cross-section reduction 
by nucleon spin fluctuations becomes large when a typi- 
cal nucleon spin-fluctuation rate is of order the ambient 
temperature T or larger. Since nucleons interact by a 
spin-dependent force, the spin-fluctuation rate is roughly 
identical with the nucleon collision rate. Therefore, the 
condition that the spin-fluctuation rate be much less than 
T ensures both that the weak interaction rates are not 
much affected by nucleon spin fluctuations and that Lan- 
dau's condition is satisfied. 

The main problem in the degenerate case is to identify 
the quantity which is to be interpreted as the relevant 
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effective spin-ffuctuation rate. Because oniy the spins of 
nucleons near the Fermi surface have a chance of evolv- 
ing in a nontrivial way, and because Landau's condition 
pertains to the quasiparticles near the Fermi surface, it 
is clear that we need to define an appropriate effective 
spin-fluctuation rate for the quasiparticles near the Fermi 
surface. 

The impact of nucleon-nucleon collisions on weak inter- 
action rates is best understood in the language of linear- 
response theory where the medium is described by the 
dynamical density and (iso) spin-density structure func- 
tions. This method allows for a straightforward calcula- 
tion of the reduction of weak interaction rates in the per- 
turbative limit where the Landau condition is fulfilled, 
and thus allows for a delineation of the physical param- 
eters where this treatment is justified. 

We will limit ourselves to the simple situation of a 
nonrelativistic, single-species medium, i.e. we will study 
nonrelativistic degenerate neutron matter. This excludes 
the important Urea processes from consideration which 
are more difficult to address because they involve two 
degenerate Fermi seas (protons and neutrons) with vastly 
different Fermi momenta. We believe that for the present 
exploratory purposes a simple toy model is best suited to 
illuminate the issues at hand. Therefore, we shall limit 
our attention to the neutral-current scattering process 
v+n — > n+is in nonrelativistic degenerate neutron matter 
in the presence of interactions which cause the neutron 
spins to fluctuate. 

In Sec. II we introduce the relevant spin-density struc- 
ture function and derive a simple sum rule which is used 
in Sec. Ill to calculate the vn scattering cross-section 
reduction by neutron spin fluctuations. In Sec. IV we 
summarize and discuss our result. 



II. DYNAMICAL STRUCTURE FUNCTIONS 

A. Definition 

In nonrelativistic neutron matter all weak interaction 
rates are determined by the dynamical density and spin- 
density structure functions. In an isotropic medium they 
are given by Q 

1 f + °° 

S p (u,k) = — / dte Mt (p(t,k)p(0,-k)}, 
nBV J_ oc 

a r+°° 

S a (u;,k) = —— dte^{a(t,k)-tr(0,-lL)), (1) 
Zn B V J_ OQ 

where tib is the baryon (here neutron) density, V the 
volume of the system, k the momentum transfer, and lj 
the energy transfer from the weak probe to the neutron 
medium. Further, p(t, k) is the spatial Fourier trans- 
form at time t of the neutron density operator p{x) = 
tp^ (x)ip(x) where ^p(x) is the neutron field operator, a 



Pauli two-spinor. Similarly, er(t, k) is the Fourier trans- 
form of the spin-density operator cr(x) = (x)Tip(x) 
with t a vector of Pauli matrices. The expectation value 
(. . .) is taken over a thermal ensemble so that detailed 
balance S p , a (ui, k) = S p ^(—uj, — k)e w / T is satisfied. Note 
that a positive lu is energy given to the medium by the 
weak probe. 

The energies of the neutrinos which interact with the 
medium are much smaller than the neutron mass or mo- 
menta so that the long-wavelength limit k — > is an 
adequate first approximation. In practice, its validity 
is questionable if neutron-neutron correlations or collec- 
tive modes are important which for the moment we shall 
assume is not the case. Then the neutrino differential 
scattering cross section is given by 

da _ Gl t\ ( 2 S p (ei - e 2 ) 2 S t7 (e 1 - e 2 ) \ 

dT 2 -^r [ Cy 2^ +6Ca 2^ J' 

(2) 

where £1,2 is the energy of the incoming and out- 
going neutrino, respectively, and S Pj<t (uj) stands for 
limk-fO Sp,o{w, k). Further, Gf is the Fermi constant, 
and Cy = — 1 and Ca ~ —1.15 are the neutral-current 
weak coupling constants for the neutron [0] . In bulk nu- 
clear matter, Ca may be suppressed somewhat. 

In a nonintcracting medium, the density and spin- 
density operators remain constant so that the dynami- 
cal structure functions are proportional to S(uj). In the 
nondegenerate case, they are S p>(T (ui) = 2ird(u>). To in- 
clude neutron-neutron anticorrelations induced by the 
Pauli exclusion principle one evaluates the expectation 
values in Eq. (|l|) by normal ordering of the neutron field 
operators, taking proper account of the anticommutation 
relations. Then one arrives at the intuitive result 

Sp ,4cj) = 2nS(cj)±Jj0f p (l-U), (3) 

where / p is the occupation number of the neutron field 
mode p. In the nondegenerate limit one may neglect the 
Pauli blocking factor (1 — / p ) so that one arrives at the 
previous result if one notes that tib = J ./p2d 3 p/(27r) 3 . 
Here, the factor 2 counts the two neutron spin degrees of 
freedom. 

Even after "turning on" interactions between the neu- 
trons, or between the neutrons and some external poten- 
tial, the density operator remains constant. The vector 
current quantity that does fluctuate in the presence of 
interactions is the neutron velocity which in the nonrela- 
tivistic limit is small. Therefore, S p (uj) remains propor- 
tional to S(uj). 

However, if the interaction involves a spin-dependent 
force as expected for neutron-neutron interactions, the 
spin-density structure function will be broadened because 
the spin of a given neutron near the Fermi surface will 
"forget" its initial orientation roughly after the collision 
time t. The width of S a (u>) roughly represents r _1 so 
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that the Landau condition corresponds to the require- 
ment that the width of So-(lo) must be much less than T. 
If this is satisfied, the neutrino scattering rates and thus 
the neutrino opacities are well approximated by the non- 
interacting result for S a (cu). Of course, it may be modi- 
fied by neutron- neutron correlations or collective modes, 
effects that were the main focus of many of the previous 
papers |f). 



is the degeneracy parameter in the nonrelativistic and 
very degenerate limit with Ep = p F /2m the nonrelativis- 
tic Fermi energy. 

Therefore, in a noninteracting degenerate medium the 
structure functions are S p ^{lo) = (3/2rj) 2tt5(lo). The 
total scattering cross section of a neutrino with energy 



ex is then a = (3/2t?) (C£ + 3C 2 A ) G| e 2 /47r. 



B. Normalization 

An important general property of the dynamical struc- 
ture functions is their normalization. If one integrates 
both sides of Eq. ([!]) over dtu, the term e ZUJt yields S(t) 
so that the time integral can be trivially done. Then the 
normalization for the spin-density case is 
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r(0,0)-<x(0,0)). 



(4) 



If one ignores spin-spin correlations, the r.h.s. is inde- 
pendent of the neutron spins' evolution. For the sake of 
argument one may imagine that this evolution is caused 
by the interaction with some external potential rather 
than by neutron-neutron collisions so that there is no 
reason to expect spin-spin correlations. 

In this case one may evaluate the r.h.s. of Eq. (|4|) as 
above and finds 
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(5) 



The occupation numbers are given by a Fermi-Dirac dis- 
tribution so that the r.h.s. is 
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e (S-/i)/T + 1 



(6) 



where E = p 2 /2m is the kinetic energy, m the neu- 
tron quasiparticle effective mass, and fi the nonrelativis- 
tic neutron chemical potential. Then Eq. (0) is 
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dpp 



(e» + l) 2 



(J) 



where z = (E — f^)/T. For very degenerate conditions 
the integrand is strongly peaked near z — (the Fermi 
surface) so that after a transformation of the integration 
variable to z one may replace p with pp and one may 
extend the lower limit of integration to —00. The integral 
can then be evaluated analytically so that altogether 



+ °° du> 3 



2tt 



2 V 



Here. 



Ep 
T 



2mT 



(8) 



(9) 



C. Perturbative Representation of Sa-(u>) 

If neutrons interact by a spin-dependent force it causes 
a nontrivial evolution of their spins and thus a nonzero 
width of 5,7 (w). Except in the neighborhood of lo = 
where multiple-scattering effects become important, 
Sa{ui) can be calculated on the basis of a bremsstrah- 
lung or medium-excitation amplitude (^fl. Because for 
small u> the result generically varies as w it is useful to 
represent it in the form 



^ ) = ?±e(u,/T)x 

■ 1* 



1 



J I T for u) < 0, 
for lo > 0. 



(10) 



The explicit distinction between positive and negative en- 
ergy transfers represents the detailed-balance condition. 
Further, s(x) is an even function which is normalized 
such that s(0) = 1. In the classical limit of hard col- 
lisions one has s(x) = 1 for all x as discussed in more 
detail in Refs. ^,^]. In general, s(x) embodies informa- 
tion about the detailed form of the interaction potential 
and about quantum corrections to the classical result. In 
the nondegenerate case, r CT has the interpretation of an 
average spin rate of change, or conversely, T^ 1 is the ap- 
proximate time for a given nucleon spin to relax, i.e. to 
forget its initial orientation. 

Explicit calculations of T a and s(x) exist for a single- 
species nuclear medium where the nucleon interaction is 
modelled by a one-pion exchange (OPE) potential [||. 
For a degenerate medium the relevant expressions can be 
extracted from Ref. 101 



r CT ,OPE = 47ra 2 r T 3 /p!, 



(11) 



where the neutron Fermi momentum is given by 115 = 
P^/3tt 2 , = (/2m/m 7r ) 2 /47r w 15 with / w 1.0 is the 
pion fine-structure constant, m is the neutron mass, and 
the pion mass has been neglected in the OPE potential. 
One also finds from Ref. 



sope(z) 



(x 2 + 4tt 2 ) |x| 



4tt 2 (1 



')' 



(12) 



which is 1 at x = while for \x\ ^ 1 it is |a;| 3 /47r 2 . 

Sigl § has derived an f-sum rule which implies that 
the integral J S a (cj) lo dio must exist and thus that s(x) 
must be a decreasing function for large x. This conclu- 
sion also pertains to the degenerate case: if the energy 
transfer to to the medium far exceeds the Fermi energy 
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Ep, a nucleoli is lifted far above the Fermi surface so 
that degeneracy effects cannot cause a modification of 
the nondegenerate result. Thus, the degenerate and non- 
degenerate s(x) must be identical for \x\ 3> Ep/T apart 
from a multiplicative factor which arises because of our 
normalization s(0) = 1. 

Explicit calculations of s{x) for various assumptions 
concerning the neutron interaction potential and for var- 
ious degrees of neutron degeneracy are left for a future 
study @. 



D. Physical Estimate of T„ 

It will turn out that the vn scattering cross-section 
reduction is primarily sensitive to the neutron spin- 
fluctuation rate r CT . Therefore, it is useful to understand 
on physical grounds its overall magnitude and its scaling 
with temperature and density. To this end we assume 
that neutrons scatter with a velocity-independent cross 
section a n which is caused by a spin-dependent force such 
that the neutron spin is flipped in a typical collision. If 
the interaction is approximated by an OPE potential, on 
dimensional grounds the cross section is estimated to be 
a n ~ a 2 /to 2 . We will assume that the scattering is ei- 
ther due to a random collection of external scattering 
centers with a density n c , or due to collisions with the 
other neutrons with a density ns. 

If the neutrons are nondegenerate they move with a 
typical thermal velocity v « (3T '/to) 1 / 2 . By assumption 
the spin-fluctuation rate is roughly equivalent to the col- 
lision rate so that T a w n c {a n v) w n c a n (3T '/to) 1 / 2 . With 
the above estimate for a n and with the other neutrons 
being the scattering centers (n c — ub) one finds that 
r CT scales as a 2 

T i/2 m - 5 /2^ This 

agrees with an explicit 
calculation which yields Ay/n for the numerical factor 

Next, we consider degenerate neutrons for which a typ- 
ical velocity is pp/m. If they interact with external scat- 
tering centers, the collision rate for neutrons near the 
Fermi surface is about n c a n pp /to. However, only the 
scattering of neutrons with an energy E within about a 
distance T from the Fermi surface is not blocked by de- 
generacy effects. This is an approximate fraction T / E-p = 
I/77 of all neutrons. Therefore, the spin- fluctuation rate 
averaged over all neutrons is T a » n c o n {pp / m)(T / Ep) . 

Finally, if the scattering is among degenerate neutrons 
we have n c — ns and a typical relative velocity pp/m. 
The average collision rate among neutrons is reduced by 
several factors of ?y _1 = T/Ep. Two such factors arise 
because each initial-state neutron must have an energy 
within about T of the Fermi surface. One further factor 
arises because one final-state neutron must also lie near 
the Fermi-surface; energy-momentum conservation then 
ensures that the other one fulfills this condition as well. 
Altogether we thus find r CT sa <j n nB(pp/m)(T/Ep) 3 . 
With ns = pp/3n 2 and a n « a 2 /m 2 we thus recover 
Eq. (O) apart from the numerical coefficient. This T a is 



the spin-fluctuation rate averaged over all neutrons. The 
spin-fluctuation rate of those neutrons which lie near the 
Fermi surface is larger by a factor n. 



III. CROSS-SECTION REDUCTION 

A. General Result 

We may now proceed to calculate the vn scattering 
cross section in the presence of spin fluctuations of the 
degenerate neutrons. To this end we begin with the total 
axial-current scattering cross section a a of a neutrino 
with energy e±. In the structure-function language it is 
the dei integral of the axial part of Eq. (||) or equivalently 
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^H( £l - W ) 2 e( ei - W ). (13) 



The problem with this expression is that it diverges if one 
uses the perturbative expression S^ rems (u>) instead of the 
full but unknown Sa-hv). Following the treatment of the 
nondegenerate case || we note, however, that Eq. jl^ ) 
can still be evaluated on the basis of S^. Ioms (ui) without 
knowledge of the detailed low-u; behavior if one includes 
the sum rule Eq. (||). 

To this end we note that for degenerate free neu- 
trons the vn scattering cross section is OA,iree = 
(3/2^) (3/47r)C^G|e 2 as stressed after Eq. ©.'There- 
fore, the interaction-induced modification 8a a = oa — 
CTA,froc is given by 



$cr A 

CA,fre< 



= -1 



du 2n S a (uj) (ei 
2^ 3 



-w) 2 e( £l -w) 



(14) 



Then we may proceed analogously to Ref . |3j and replace 
— 1 with an integral over the structure function by virtue 
of the sum rule Eq. (ph, 



Scr A 

OA, fre< 
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du 2rj Sa-(u>) 
2^ J! 



(e!-^) 2 ©^-^) 



(15) 



For small u> the integrand varies effectively as ui 2 S a (ui) 
because the term linear in oj switches sign at the origin. 
Therefore, to lowest order we may substitute ^(w) — > 
S^ Tems (oj), provided we interpret the remaining integral 
by its principal part. 

This result becomes more transparent if we consider 
the reduction of an average cross section rather than 
one for a fixed initial-state neutrino energy. To this 
end we use a Maxwell-Boltzmann distribution of neu- 
trino energies at the same temperature T which char- 
acterizes the ambient neutron medium. The thermally 
averaged free cross section is found to be (aA,free) = 
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dxx 2 G(x)s(x). 



(18) 



(3/277) (9/ir)ClGlT 2 . Because Eq. @ is fully analo- 
gous to the corresponding result of Ref. §| apart from 
an overall factor 277 /3 we may conclude without further 
calculations that 

5{a A ) 2n f 00 dx ~ 

7 \=-^r r Sf W G W' ( 16 ) 

\0"A,free/ o Jo ^ 7r 

where x = u>/T, 

G(x) = 1 - (1 + x + \ x 2 ) e- x = \x 2 + 0(x 3 ), (17) 

and S„(x) = TS a (xT). 

As in Ref. || the x 2 behavior of G{x) at small x allows 
us to replace S a (x) to lowest order with the perturbative 
S^ iems (x). Therefore, with the representation Eq. jic| ) 
and with 7^ = T a /T we find for the cross-section reduc- 
tion 

(f A, free) 3 27T 

The integral expression is 5/6 for the classical approxima- 
tion six) — 1. In general, the integral will be a numerical 
expression of order unity. Its precise value for a variety 
of assumptions concerning the cause for the neutron spin 
fluctuations will be studied elsewhere JlO| . 

Equation ( |l8| ) shows that the expansion parameter 
which defines the perturbative regime is 2rj^ a /?>, as op- 
posed to the nondegenerate case where it was found to be 
7 ff . In both cases 7 CT is defined to be the spin-fluctuation 
rate averaged over all neutrons of the medium. How- 
ever, in the degenerate case only the neutrons near the 
Fermi surface participate in collisions; it is their spin- 
fluctuation rate which reduces the vn scattering cross 
section. The quantity 2nj a /3 corrects for this effect. It 
is to be interpreted as an effective spin-fluctuation rate 
for the neutrons near the Fermi surface, in agreement 
with our estimates of Sec. II. D. 

We conclude that a "naive" perturbative calculation 
of neutrino interaction rates in a degenerate neutron 
medium is possible if nY a <C T while significant correc- 
tion arise if nT a > T, This latter case corresponds to 
a situation where the collision rate of neutrons near the 
Fermi surface is not small relative to T, in violation of 
Landau's condition for the applicability of a Fermi-liquid 
treatment. 

In the nondegenerate case it was reasonable to extra- 
polate the behavior of the cross section (a A ) into the 
nonperturbative regime by virtue of an explicit ansatz 
for the low-w behavior of S a (uj) which incorporated the 
equivalent of the sum rule Eq. (ph. In the present case 
such an extrapolation is far more problematic because 
the derivation of the sum rule itself was based on the 
assumption that neutrons can be treated as quasiparticles 
which follow a thermal Fermi-Dirac distribution. In the 
nonperturbative regime this assumption is not justified 
so that in the present case the sum rule has a weaker 
standing than it did in the nondegenerate case where we 
did not need to invoke the anticommutation relations for 
the nuclcon fields. 



B. Numerical Result for OPE Potential 

If neutron-neutron collisions are the primary cause for 
neutron spin fluctuations, and if one models the inter- 
action by an OPE potential, we may use Eq. (11) to 
estimate T a . In this case we find 



2?7 r^ppE 
3 T 



4tt 



T 



T 



T a; m =L0 °MeV' 



(19) 



where we have used the vacuum neutron mass for the 
numerical estimate. This result does not depend on the 
density (or Fermi momentum) which fortuitously cancels 
as explained by the physical arguments in Sec. II. D. If 
the neutron spin fluctuations were caused by the interac- 
tion with a distribution of external scattering centers, T a 
would depend on their density as well as on the neutron 
Fermi momentum. 

Of course, if neutron-neutron interactions are the pri- 
mary cause for neutron spin fluctuations one would also 
expect significant spin-spin correlations which we have ig- 
nored. However, in order to study spin-spin correlations 
in the framework of a Fermi liquid theory one would need 
to assume that Landau's condition is fulfilled which is not 
the case in any obvious sense when nr a > T. Therefore, 
we believe that the usual calculations of neutrino opaci- 
ties in hot degenerate neutron-star matter are applicable 
only for T < 1 MeV. 



IV. SUMMARY 

We have derived an expression for the vn scatter- 
ing cross-section reduction in degenerate neutron mat- 
ter caused by neutron spin fluctuations. We have used 
the linear-response theory approach of Ref. || , but un- 
doubtedly one would reach the same result by the direct 
perturbative method of Ref. [Q . 

In a neutron star, these spin fluctuations will be caused 
by a spin-dependent nn interaction potential. Therefore, 
in general spin-spin correlations will also be important 
which may cause further modifications of the scattering 
cross section. Many of the previous papers which deal 
with weak interaction rates in neutron stars were ded- 
icated to an analysis of these latter effects. We stress, 
however, that these calculations were based on the as- 
sumption that Landau's condition is satisfied which is 
roughly equivalent to the requirement that the autocor- 
relation function of a single nucleon spin near the Fermi 
surface is narrow on a scale set by the ambient tempera- 
ture T. 

If neutron-neutron interactions are modelled by a one- 
pion exchange potential we estimate that the usual per- 
turbative calculations are justified for T < 1 MeV, a tem- 
perature which is reached very quickly in a neutron star 
after formation. Therefore, the long-term cooling history 
remains unaffected. Of course, a calculation of the long- 
term cooling history does not require knowledge of the 
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neutrino opacity anyway as at late times neutrinos are 
no longer trapped. Roughly speaking, then, the neutrino 
opacities matter only for a short period after formation of 
a neutron star. However, precisely for this period the vn 
scattering rate cannot be calculated by straightforward 
perturbative techniques on the basis of first principles. 
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